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Quantum field theory is the application of quantum physics to fields. It provides a theoretical 
framework widely used in particle physics and condensed matter physics. One of the most distinct 
features of quantum physics with respect to classical physics is entanglement or the existence of 
strong correlations between subsystems that can even be spacelike separated. In quantum fields, 
observables restricted to a region of space define a subsystem. While there are proofs on the 
existence of local observables that would allow a violation of Bells inequalities in the vacuum states of 
quantum fields as well as some explicit but technically demanding schemes requiring an extreme fine- 
tuning of the interaction between the fields and detectors, an experimentally accessible entanglement 
witness for quantum fields is still missing. Here we introduce smeared field operators which allow 
reducing the vacuum to a system of two effective bosonic modes. The introduction of such collective 
observables is motivated by the fact that no physical probe has access to fields in single spatial 
(mathematical) points but rather smeared over finite volumes. We first give explicit collective 
observables whose correlations reveal vacuum entanglement in the Klein-Gordon field. We then 
show that the critical distance between the two regions of space above which two effective bosonic 
modes become separable is of the order of the Compton wavelength of the particle corresponding to 
the massive Klein-Gordon field. 

PACS numbers: 03.70.+k, 03.67.Mn, 03.67.-a 



I. INTRODUCTION 

Entangled states of composed quantum systems are a subject of particular interest since they 
manifest genuinely quantum mechanical properties — they yield correlations between observables 
measured on the subsystems that cannot be explained by any local realistic model 1] . Since entan- 
glement is the primary resource that allows quantum communication and computation protocols to 
outperform classical ones Q , its investigation reaches far beyond fundamental concepts of quantum 
physics. It also has a central role in some macroscopic properties of solids - internal energy, heat 
capacity or magnetic susceptibility can reveal the existence of entanglement within solids in the 
thermodynamical limit (34a] - Entanglement is further related to superfluidity, the Meissner effect 
and flux quantization in superconductors as well as long range order correlations in Bosc-Einstein 
condensates (BEC) (fj. 

A natural framework for considering systems composed of parts, which are associated with dis- 
connected regions of space, is quantum field theory with its causal and local structure, where one 
may treat fields supported on spacelike separated regions as subsystems. Here we will consider 
the entanglement of the relativistic vacuum state. This is not only thought to be connected with 
black holes thermodynamics Q or the holographic principle @ , but also to manifest itself in the 
Bekenstein-Hawking black hole radiation [t| and Unruh acceleration effects [l(| • 

There are a number of studies proving entanglement between spatial regions in the ground state 
of relativistic quantum field theory. The central result underlying all the "proofs in principle" is 
the Reeh-Schlieder theorem [1JJ, formulated in the language of algebraic quantum field theory. Let 
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A be a space-time region and A(A) the local algebra of all the operators with support in A Q. The 
theorem states that, for any state of finite energy l^) (in particular, the vacuum) the subspace 
A(A)\ip) is dense in the entire Hilbert space, which means that for an arbitrary state \tp'), we 
can find a sequence of local operations IT ra such that lim Il _ i . 00 H n \i/}) = \ip'). In particular, we can 
reconstruct an arbitrary state \ipB) in region B by applying such operations in A and then tracing 
outside B. The vacuum is thus an entangled state as it allows for a remote state preparation jl 2j j . 

As a consequence of the Reeh-Schlieder theorem and the positivity of a partially transposed 
separable density matrix, it is possible to prove that causally separated local regions are entangled 
in the vacuum state [lj]. There exist operators A, B in algebras A{A), B(B) associated with 
causally separated regions A, B that allow for constructing an entanglement witness (an operator 
whose mean value is larger or equal to zero for all separable states) , which is violated (i.e. negative) 
in the vacuum reduced to these two regions. It is possible to choose operators A, B as annihilation 
or creation operators of some local bosonic modes but the explicit construction is not known. 

It is also possible to prove the existence of local, bipartite observables that allow violation of Bell 
inequalities in the vacuum state (l4| . but, again, the specific form of these observables is not known. 
Another result is that it is possible to locally couple two qubits (detectors) to the field in such a way 
that, after a finite time, the reduced state of such a pair of two-level systems has negative partial 
transpose [Hj]. In this last case, one needs an explicit expression for the space-time dependence 
of the detector-field couplings. In order to prove entanglement for arbitrary separation of the 
detectors, these couplings need to be of a very specific, fine-tuned, form, involving superoscillating 
functions that require switching the sign of the interaction between the field and probes during the 
experiment. This seems technically extremely demanding. Along different lines (without referring 
to observables) in 16] entanglement between two separated segments of one-dimensional free Klein- 
Gordon field in the vacuum state is quantified by the logarithmic negativity which is investigated 
numerically and shown to be finite for both critical and noncritical field limit. 

To which extent can the vacuum of a quantum field be operationally accessed and used as an 
entanglement resource? By integrating over scalar Klein-Gordon field operators with compactly 
supported real functions - i.e. detection profiles localized in two regions of space - we define col- 
lective field operators (and collective conjugate momenta). These weighted averages of operators 
allow reducing the vacuum to a system of two effective bosonic modes. To such a reduced state 
we apply the entanglement measure for continuous variables systems based on Simon's criterion 
for separability [l7j . This approach has several advantages. Since separability criteria for infinite- 
mode states are unknown, we need to reduce the vacuum to a more comprehensible system. Our 
approach enables to quantify entanglement present in the resulting two-mode system, which is still 
infinite dimensional, unlike in (l5j . where entanglement is first transferred to two qubits. Physi- 
cal probes have finite spatial resolution, so introducing collective observables is a reasonable first 
approximation towards a more realistic treatment of the problem. Finally, the Klein-Gordon field 
is the continuum limit of an infinite linear harmonic chain, and within this "collective approach" 
entanglement between blocks of oscillators in the ground state of the chain was found and quantified 
in [l8[ . It is therefore interesting in itself to understand the relation between the field and the chain 
from this particular perspective. 

In this paper we prove the existence of a critical distance between two regions of space above which 
two effective bosonic modes associated with the regions become separable. From the numerical 
analysis of the linear harmonic chain this critical distance is estimated to be of the order of the 
Compton wavelength of the particle corresponding to the massive one-dimenisonal Klein-Gordon 
field (the continuum limit of the linear harmonic chain). We also give an explicit example of the 
possible profiles that allow for a construction of entangled modes. Numerical results obtained for 



1 More precisely, A(A) is defined as the algebra generated by operators of the form f A dxdtf(x,t)Q(x,t), 
J A dxdt g(x,t)n(x,t) with test functions /, g. 
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this exemplary functions are presented and discussed. 

The structure of the paper is as follows: In Sec. [IT] the entanglement criterion and measure, which 
are further applied, are introduced and commented. Section|nT]is devoted to the collective operators 
in the relativistic scalar quantum field theory framework. We give their definition in terms of the 
detection profiles and discuss the constraints on the latter. Next, we derive one of the main results: 
the proof of the existence of a critical distance above which the modes defined within the collective 
operators approach, become separable. The section is closed with an explicit construction of the 
interaction needed to read out our observables. Section IIVI deals with the linear harmonic chain. 
First, the relation between harmonic chain and continuous field is briefly reviewed. Second, on 
the basis of numerical analysis, the optimal profiles (maximizing the entanglement witness) in the 
discrete case are described and the critical distance for them is obtained. For the optimal profiles 
in the continuum limit this distance is then estimated to be of the order of the particles' Compton 
wavelength. In Sec. [V] we give an explicit example of profiles that allow for defining entangled 
modes and present the numerical results obtained for this special case. We close the paper with 
final remarks and conclusion in Sec. IVI1 

II. SEPARABILITY CRITERIA 

In this paper we use a particular form of the separability criterion derived by Simon [l7j , necessary 
and sufficient for two-mode Gaussian states. Following the original notation, we introduce a vector 
of phase space operators for the two modes system: 

|= (Qa,Pa,Qb,Pb) T 

where T denotes transposition. Canonical commutation relations (CCR) in natural units, i.e. 
h = c = 1 , can be concisely written in a matrix form 

[&> = i^ij > 
where we use the two-mode symplectic matrix 

°»(.°. ;)©(-! iy 

Defining the variance matrix Vy := | — (£i),£j — (£j)}^, where {A, B} := A ■ B + B ■ A, we 
obtain a compact statement of the Heisenberg uncertainty relations: 

v+ l -n>o. (1) 

Every physical state has to satisfy this inequality. For separable states it must hold also after 
partial transposition 0. The effect of partial transposition at the level of variance matrix elements 
is only (PaPb) — > — (PaPb), provided that the system also satisfies 

Al. (Qi) = (Pi) = , ie {A, B} 
A2. ^{{Qi,Pj}) = 0, i,j G {A,B}, 



2 For relativistic quantum field theory partial transposition on the level of field operators is a partial CPT (completely 
positive and trace preserving)-map combined with reflection [13| |. In our case it reduces to partial time reversal, 
exactly as in nonrelativistic continuous variables systems [17| . 
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which is always the case for us (one can always find a local symplectic transformation that enforces 
Al and A2 [H Gl). Thus > a11 P h y sical 

separable states fulfilling Al, A2, satisfy the following 

inequality (Simon's criterion): 

I (Qa)(P1) - (Q 2 b)(P 2 b) - 2\(QaQb)(PaPb)\ + 

(2) 

+4 ((QIXQI) - (QaQb) 2 ) ((P 2 a )(P 2 b ) (PaPb) 2 ) > o . 
In the case that the system additionally satisfies 

A3. (Q\) = (Q%) , (P 2 A ) = (Pi) , 
we can considerably simplify the separability criterion by first factorizing it 

{{- - ((Q 2 a) - \(QaQb)\) ((P 2 a ) - \(PaPb)\)}x 

x{i - ({QD + KQaQb)\) ((P 2 b ) + \(PaPb)\)}>0. 

Notice that the second factor is always non positive. If it equals zero, both correlations (QaQb) 
and (P a Pb) must vanish, but then also the first factor equals zero. So finally, the above inequality 
is equivalent to 

({QD - \(QaQb)\)-((P 2 a ) \{PaPb)\) >\. (3) 

Condition A3 is not satisfied for general profiles, but can be justified by physical assumptions (e.g., 
using the same detectors in regions A and B). In Sec. HVI this assumption will also be supported by 
numerical results. We also have some evidence that profiles that minimize the left-hand side of © 
satisfy A3 (see SecEE}. 

To quantify entanglement, we will use the degree of entanglement given by 

e=l- ±({Q\) - \{Q A Q B )\)-({Pl) - \{P A P B )\) . (4) 

For Gaussian states e > iff the state is entangled. In this case £ is a monotonically increasing 
function of the negativity N (the absolute sum of the negative eigenvalues of the partially transposed 
density matrix): e = j^rr ■ The negativity is based on the Peres-Horodecki criterion [2(| [21J and 

was shown to be an entanglement monotone [22], [23| . 

In the literature there are plenty of other criteria which, as well as the one we use, are necessary 
and sufficient for Gaussian states and also are phrased with second order correlations between 
phase space operators (e.g. [l9j], (24|, [1^|). It is then obvious that for Gaussian states they all 
yield the same results. However, for non-Gaussian states those criteria are only sufficient, so we 
should ask which of these is in general the strongest. In other words: are there any non-Gaussian 
entangled states such that one criterion out of those mentioned above would detect it, whereas some 
other would fail? The answer is negative: either all of them will be satisfied or all violated. From 
a physical point of view this is clear, simply because we cannot discriminate Gaussian from non- 
Gaussian states on the basis of their variance matrices (second order correlations) . Mathematically, 
it can be shown that all those criteria are given by the same inequality up to a local linear canonical 
transformation of modes. This means that the choice of any of these criteria results in the same 
set of states detected as entangled. In any case, as we are going to see, in this work we will only be 
concerned with Gaussian states, so that ([2]) (and (|3]), when A3 is met) will always provide necessary 
and sufficient conditions for separability. 
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III. COLLECTIVE OPERATORS 

We consider a massive Klein-Gordon (KG) field. The field (j>(x, t) and conjugate momentum 
Ti(x,t) satisfy the equal-time canonical commutation relations, [26j 

[tt(x,0, 4>(y,t)} = i5 3 (x-y) 

[$(x,t)J(y,t)]=0 (5) 
[*{£,*), 7r(y,t)}=0 

They can be expanded in terms of creation and annihilation operators, d^ and at of normal 
modes: 

i /- + 00 

4>{x,t) = —-^r / d 3 k ./-!- (aj:e tkS - mkt + tf e -^+™*t 



t) = ^ jT~ d»* ^ ( %e ^-* - ate-*"-*) 



+ 00 



where ccjfe = \/fc 2 + to 2 . 

The vacuum state is defined by the property 

a^lO) = Vfc . (6) 

We study the possibility to detect entanglement in the vacuum state when the allowed measure- 
ments are constrained to field operators smeared over two bounded space-time regions [collective 
field observables). Our motivation is twofold: first, a field operator in a single space-time point is 
not a physical observable but a purely mathematical concept. Physical probes always have finite 
spatial resolution, so introducing collective observables is a reasonable first approximation towards 
a more realistic treatment of the problem. Second, since separability criteria for infinite-mode states 
are unknown, we want to reduce the vacuum to a system consisting of only two bosonic modes. 
In a general case, the smearing is given by two different real functions gA^x), gB^x) with compact 
supports and collective operators in the regions A, B are defined as follows: 

/+oo 
d 3 xg A / B {x- x A/B )(j>(x,i) 

,+S , (7) 
P{xa/b) = Pa/b ■= / d 3 xg A/B (x~ x A / B )n(x,t) . 

J — oo 

We consider collective operators that satisfy CCR 

[Qa, PA]=i, 

[Qa, Pb}=0, (8) 
etc ... 

which hold if the profiles satisfy orthonormality conditions (to be discussed in the next paragraph). 
The fact that collective phase space operators satisfy CCR guarantees that indeed we deal with 
two distinct bosonic modes, so it is meaningful to treat them as subsystems and speak about 
entanglement or separability of their joint statc0. 



3 More precisely, the collective operators on A and B generate two commuting subalgebras, that in turn induce two 
subsystems in the Hilbert space [27|, each of which is isomorphic to the space of a one-dimensional particle. 
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Completing the set {gA(x — %a), 9b(x — %b)} to an orthonormal basis in L 2 (R 3 ), can be 
extended to a linear canonical transformation of modes, two of which coincide with our collective 
ones. Tracing the global ground state over all but these particular two modes gives the final 
state of the two subsystems. Since the vacuum of a scalar quantum field is Gaussian and both 
transformations preserve this property, the final reduced state is Gaussian as well. This observation 
is very important as it means that the criterion §5§ is in our case necessary and sufficient (if condition 
A3 is satisfied) and ((4]) indeed gives the degree of entanglement between the two bosonic modes. It 
is clear that the expectation values of all the possible products and combinations of the collective 
operators (0 in the global vacuum (0 coincide with their values calculated with respect to the 
reduced state. 

The idea of restricting the possible measurements to a pair of collective modes is an extension 
of [lH (where the chain of harmonic oscillators is considered) to the framework of scalar quantum 
field theory (QFT) with general profiles. Our main goal is to explicitly find two profiles gA/B{%) 
such that the effective modes ([7]) are entangled. 

Conditions on the profiles. Before proceeding to prove our main results, we need explicit 
expressions for the constraints that the detection profiles gA/B(%) have to satisfy. It is useful to 
express them in terms of the Fourier transform 

i r +oc 

We require that the collective operators (0 satisfy CCR (0. All the relations involving only col- 
lective position or only momentum operators are automatically satisfied due to (0 . The remaining 
ones lead to orthonormalization of {gA(x — $a), gB(x — xb)} in L 2 (M. 3 ). With i,j — A,B, we have 

/+oo 
d 3 ke- ik ^-^ gi (k) gj (-k) = 5 l3 (9) 
-co 

Further, all correlations in (0 should be finite. From the Cauchy-Schawrz inequality \(QaQb)\ < 
\j{Q A ){Q%) and \ {P a Pb) \ < \J {P\){P B )- So > ( witn assumption A3) it is only necessary that (Q 2 A ) 

and (Pi) are finite. This indeed holds, if the Fourier transforms of the smearing functions decay 
fast enoug ifi 

\<)A/B(k)\ < — r for | & | — > oo , A > — — — in d space dimensions. (10) 
\k\ 2 

Finally, we demand A1-A3. Note that Al, A2 are satisfied for every profile, while A3 is equivalent 
to the additional conditions 




4 This condition was not satisfied in a paragraph devoted to scalar quantum field in [ljjj . However, this fact does 
not affect the validity of the results obtained there for the linear harmonic chain 
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All correlations appearing in the criterion ([3j) in terms of the Fourier transforms of the profiles 
read = A, B) 

/+°° i _ 

d 3 k— e- ik &-^ gi (k) gj (-k), 
-oo lvJ k /-ii-v 

(PiPj) = / d 3 k^e- ik ^-^ 9i (k) gj (-k). 

J-oo * 

We will further denote D := x*b — x*a- This parameter, appearing in (QaQb) and (PaPb), is a 
measure of the distance (separation) between the subsystems. 

Large separations limit. It is natural to ask whether there exist profiles satisfying all the 
constraints given above and defining entangled modes for arbitrary separations. In this paragraph 
we will show that this is not possible. Below we prove that for every pair of allowed profiles there 
exists a finite critical distance D cr i t such that corresponding collective operators are separable at 
distances larger than the critical one. More precisely: for every pair of orthonormal functions 
9a/b(x), satisfying A3 and (fTU)) [giving finite correlations in @] there exists D crit < oo such that 
the modes defined by g A (x — xa), 3b (a; — xb) are separable for every \D\ > D crit . We will give the 
proof in three space dimensions but it remains valid in arbitrary finite dimensions. 

First, notice that the integral form of the Cauchy-Schwarz inequality (i.e. Holder's inequality) 
together with condition (|TU| . not only guarantees the finiteness of all the correlations (jlip but, also 
asserts that their integrands are functions from L 1 (R 3 ) (space of functions, which absolute value 

is integrable). In particular, it enables to prove that both gA ^ gB ^" k ^> an d Wkg A (k)gB(—k) are in 

From Riemann-Lebesgue lemma [28| (it says that the Fourier transform of an L 1 function vanishes 
at infinity) it now immediately follows that 

lim (Q a Qb)= lim / d s k — e~ ikD g A (k)g B (-k) = 0, 

|j3|->oo |5|->oo J-oo iw k 



d 3 k^e- ikD g A {k)g B (-k) = 0. 

In consequence, for \D\ — > oo the left-hand side of the criterion reduces to the product (Q\) (Pa)- 
Because of CCR imposed on the collective operators (JSJ, the Heisenberg uncertainty relation guar- 
antees that this product is always greater or equal 4 , so in the limit of infinite separation the modes 
become separable. However, we obtain much stronger result by utilizing Eq. Let us write the 

product (Qa)(Pa) as a double integral symmetrized over the integration variables 

(Qa)(p1) = \ [[ +0 ° d*kd 3 q 1 -(^ + ^) \g A (k) g A (q)\ 2 >\. 

4 J J-oo 1 \ W k rU q J 4 

The last inequality is a direct consequence of the normalization of the profiles [conditions ([5])]. the 
fact that + > 2 for all fc, q e M 3 and that the latter saturates only on the hyperplane 

|fc| = \q\, which has zero Lebesgue measure. 

Summarizing, we have shown that lim^^^ ((Q A ) - \(QaQb)\)-({P%) ~ \ (PaPb)\) > \, which 
is equivalent to 

3D crit < oo such that ([Q 2 A ) - \(Q aQ b)\) ■ {(P 2 A ) - \(P A Pb)\) > \ V|D| > D crit . (12) 
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This is, however, exactly the separability condition ([3]). So (fT2j) states that, given a pair of or- 
thonormal functions gA(x), gs{x) which satisfy (fTCTI) and A3, there exists D cr i t < oo such that once 
\xa — xb\ = \D\ > D cr u, the modes defined by gA(x — xa), 9b(x — xb) are separable. 

It is interesting, although not very relevant from the physical point of view, that this result is 
not restricted to functions with compact support, as this property has no specific role in the proof 
(only orthonormality is needed). For example, it holds for any pair of orthonormal test functions 
(element of Schwartz space: rapidly decreasing functions) satisfying A3. 



Measurement of the collective observables. A significant difference between our and 
other related works lies in the treatment of the detector. We define it only in terms of the observables 
that it measures. In this paragraph we propose a model by which a system interacting with the field 
can be used to implement the desired measurements of the collective observables. We will follow 
the method discussed in [29| . In order to measure each of the operators we thus need a different 
interaction. In general, however, to either of the effective modes A, B we couple another bosonic 
mode - a detector. Then, under the coupling of the suitable degrees of freedom of the subsystem 
and the detector, measurements made on the latter reveal the value of the corresponding observable. 

Let us consider a universal situation when observable W is measured on mode i, with i = A, B. 
The Hilbert space Hdi of the joint system comprised of the subsystem i (with Hilbert space denoted 
by Hi) and its detector (with Hilbert space Hd) is a tensor product Hdi = Hi® Hd- We introduce 
the phase space observables for the detector Qd, Pd, which satisfy canonical commutation relations 
[Qd, Pd] = i- In other words, we consider the detector's state space to be isomorphic to the Hilbert 
space of a one-dimensional particle. An interaction Hamiltonian, which allows a measurement of 
the operator W, takes the form 

Hf := aWPd , (13) 

where a is a time independent coupling constant for this particular interaction. As bases of the 
Hilbert spaces Hi, Hd we take the sets of eigenvectors of W, Qd, namely {\w)} w ^s.,{\qd)} qd eM. such 
that W\w) — w\w) and Qd\ld) — qd\qd)- (We consider here observable W with continuous spectrum 
because we are primary interested in measuring the collective observables, however the framework 
considered below applies also to operators with discrete spectrum.) The state of the subsystem i 
is mixed. In the chosen basis it can be written as 

-j-oo 

dw dzp-i(w, z)\w) (z\ (14) 

-oo 

The initial state of the detector is, in an idealized situation, a pure eigenstate \qd) of Qd- Therefore, 
prior to the interaction, the joint state po of the two systems is 

/• r + oo 

Po= dwdz pi(w,z)\w)(z\®\q d )(q d \. (15) 

J J — oo 

If the time scale of the measurement process is much smaller then that of the free evolution of the 
field and the detector, the time evolution of the density matrix po is given by the Hamiltonian (1131) : 

p t =e-^% e^ Wt . (16) 

Making use of the fact that the momentum operator is a generator of spatial translations (see e.g. 
3f|) we notice that 

e- lH ' *|w) <g> \q d ) = \w) <E> \q d + atw) , (17) 
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so finally 



1 \ / , , / W- q d z- q d 



at / ./-co V at at 



q d \ z- q d 



at / \ at 



®\w)(z\. (18) 



Performing measurements on the detector's degrees of freedom, we can reconstruct the values of 
the observable W in the state of subsystem i given by (TlH) . It is straightforward to derive that the 
probability amplitude for the detector to be at time t in some eigenstate \a) of Q d is proportional to 
Pi ( a ~at d ; ^f^)j which is in turn the probability amplitude of subsystem i being in state | a ^ d ). A 
more realistic treatment would involve assuming for the initial state of the detector not an eigenstate 
of Q d but rather a superposition J_°° dq d f(q d )\q d ) [e.g. for a coherent initial state f(q) would be a 
Gaussian packet]. In such a case the amplitude for the detector to be at time t in the state \a) is 
proportional to dq d Pl (^M, ^) \f(q d )\ 2 . 

In the above example the considered subsystem was coupled to the detector's degree of freedom 
which corresponds to the operator P d . This choice is of course arbitrary, i.e. equivalently well the 
other degree of freedom, corresponding to Q d , may be utilized to perform the measurement. In such 
a case the interaction takes form Hj W := f3WQ d , where f3 is again a time independent coupling 
constant. As a basis of the Hilbert space H d we take the set of eigenvectors of P d and, as the initial 
state of the detector, we consider the eigenstate of P d . Following all the previous steps with these 
changes in mind, we obtain that measuring this detector's state in the momentum basis, again, 
enables to reconstruct the value of the observable W on the state of subsystem i. 

Although we couple each mode with an (effective) one-dimensional particle, different degrees of 
freedom of the measured modes are involved in measurements of conjugate collective observables. 
So, in the outlined scheme the interactions are relatively simple but the detectors serve solely as 
devices to reconstruct the values of collective operators in the two-mode state considered. On the 
other hand, if the vacuum entanglement is actually transferred to the detectors (as in |15]), by 
the price of a very fine tuned, time dependent interaction, it is possible to detect entanglement for 
arbitrary separations^. 



IV. RESULTS FROM DISCRETE SYSTEMS 



The result of the previous section says that no matter what are the shapes of the detection 
profiles, collective operators too distant from each other must be separable, but it does not give 
any indication about the possibility of finding entanglement by measurements in regions sufficiently 
close together. It is not trivial to tackle this problem directly, as the possible detection profiles form 
an infinite dimensional space and many of them will still define separable modes; it is therefore 
instructive to study a discretized version of the system, where numerical analysis can be performed. 
The numerics will give important insights into the shape of the profiles that can show entanglement 
and, furthermore, it will give evidence of the existence of a critical distance independent of the 
specific profile, thus strengthening the results of Sec. IIII1 and provide an estimation for it. Here 
only the one-dimensional case will be considered. 

Continuum limit. A one-dimensional bosonic field can be formally obtained as the continuum 
limit of a chain of coupled harmonic oscillators. This is especially useful for us, as the discrete 
system can be better analyzed; in particular, the freedom in the choice of profiles reduces to an 
optimization problem of functions of a finite number of degrees of freedom. 



5 If the detector model is described by a natural interaction, it is not possible to detect entanglement with its use 
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Let us briefly review the relation between the harmonic chain and the continuous held (see, e.g., 
[12). The KG Hamiltonian 

Hkg = \ J dx (m 2 0(x) 2 + tt(x) 2 + (V0(x)) 2 ) (19) 
can be written as the limit for Ax — > of 

H dls = ^Ax£ (m 3 % + n 2 + -L(& - h-tf 

3 ^ 

= ( Aa ™ 2 $ + Ax ^ + ^.fa " ^-i) 2 ) . ( 2 °) 

where the discretized held operators are dehned as <pj := (j)(jAx), ttj := 7r(j'Ax), with j integer and 
Ax being the spacing between successive points. 

The expression (j2"0")) can be put in correspondence with a chain of N harmonically coupled 
oscillators, with conjugate observables satisfying [q~j,pk] = idjk and Hamiltonian 

H=\Y. [m^ 2 + | + Mn 2 (I, - |,_0 j , (21) 

where M is the mass of each individual oscillator, ui its proper frequency and f2 the coupling 
frequency. Periodic boundary conditions go = Qn are assumed. 
If we set (f20|) equal to (|2Tj) , we obtain 



from which we derive 



Axm 2 (j> 2 = Muj 2 q 2 , 

-2 

AxTT 2 = % , 



(22) 



mAx=% 



In order to define a correct continuum limit, the scaling of the parameters in (f2~Tj) with Ax must 
obey (|22[) with m fixed. 

In order to simplify the analysis, we can rewrite pip in the following form: 



E N 

^ = yE(^f ft w-i). ( 23 ) 

where Eq = \/2il 2 + uj 2 , a = 2£1 2 /(2J7 2 + lu 2 ) and we introduced the dimensionless variables 
q~j = Cq~j , pj = pj j C, with C = \J Mu>{\ + 2il 2 /w 2 ) 1 / 2 . In this form, the system is characterized by 
a single dimensionless parameter, the coupling constant a, which, by construction, is constrained 
to values < a < 1. In this case, the continuum limit is obtained by setting 

' (24) 



1 + ±Ax 2 m 2 
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and letting Ax — > 0, with m constant. 

If we want to describe a region of size L using ([23]) as a discrete version of a Klein-Gordon field, 
we have to consider in the chain a block with a number of sites 

L I a 

n = — = Litla / — — . 

Aa; V 2 - 2a 

As n has to be an integer, for some values of a, L and m the expression above is not well defined. 
It can therefore be more convenient, especially for carrying out numerical computations, to express 
a as a function of n and the physical length: 




This relation fixes the physical size of a region. By increasing n and having the coupling constant 
scaling as in ([25]), one approaches in the limit a region of size L of a KG field with mass m. It is 
worth stressing that it is not the number of points n that determines the size of a region, rather, this 
needs to be fixed through the relation (|25|) . Increasing n only provides a more refined description 
of the system. 



Collective entanglement and optimal profiles. The discrete version of the collective op- 
erators ((7]) can be defined for a block A of n sites in a chain: 



Qa ■= E"=i fjQj+l 
Pa '■= E?=i fjPj+i > 



(26) 



where I + 1 is the position of the first site of A and fj is the detection profile, determining how 
much each site in the chain contributes to the collective observables (note that the indices of the 
profiles always run in the range {1, . . . , n}, regardless of the position of the block in the chain). If 
the profile satisfies the normalization condition 

n 

E/' = 1 > ( 27 ) 



Qa,Pa 



i and the subsystem 



then the collective operators have canonical commutation relations 

they define is a bosonic mode. 

Two detectors placed in two regions A and B, described by detection profiles fj and gj respec- 
tively, effectively detect two bosonic modes. If the global state is Gaussian (as is the case for the 
vacuum state) then also the reduced state over the two modes is so, as was discussed in Sec. IIIII 
we can therefore apply Simon's criterion @ to establish whether the two modes are entangled 
or not. If e is the corresponding entanglement measure, e > meaning entanglement and e < 
separability, we can ask for which profiles fj, gj e is maximized when all the parameters are fixed. 
If two profiles exist such that e max > 0, then we can conclude that it is in principle possible to see 
entanglement, while e max < proves separability for all possible collective operators. 

Let us study systematically the case where the two regions A and B have the same length L. 
We want to know if, for a given separation D between A and B, it is possible in principle to find 
entanglement and what are the shapes of the profiles that maximize the entanglement measure. 
We will set from now on the mass of the field to m = 1, so that all the lengths will be expressed in 
units of Compton wavelength A c = -j-. Note that m is the only dimensional parameter in (I19[) . so 
it defines the natural scale of the system. Using the discretized field, we can approach the problem 
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numerically: we first fix the size L of the regions, the number of sites (i.e. oscillators) in each block 
n and the number of sites d separating the two blocks. The physical separation D is then given b}0 

d , s 

D = -L (28) 
n 

(see Fig. [1]) and a is determined by the relation (|25|) . 




FIG. 1. Two blocks of a harmonic chain A and B. In this example, each block consists of n = 5 oscillators 
and the blocks are separated by d = 3 oscillators. The physical length L of the blocks and their separation 
D are related by D = —L 



We consider a system of infinite total length, which means — ^> 1. We verified numerically that 
this limit is already well approximated for N = 10(2rt + d). With all the parameters fixed, the 
entanglement measure £ is a function of the In real variables {fj, 5j}™ =1 and we can then find 
the numerical extremum with both / and g subject to the constraint (l27l) . The first numerical 
evidence is that the optimal profiles are always mirroring, that is to say, e is always maximized 
by functions satisfying gj = f n +i-j] this allows us to reduce the problem to a maximization over 
n variables. Furthermore, this symmetry ensures that (Q\) = (0%), (Pa) — (Pb) (condition A3 
above), so that we can use the simplified entanglement measure (j4]). Notice, that the obvious choice 
of a rectangular profile (a "top-hat" function) is far from optimal. Moreover, such a profile would 
not work in the continuum limit, as it gives diverging correlations. 

We proceed in the following way: first we fix the physical region size L, then, for a given value 
of d, we look for the critical block size n cr n(d) such that the blocks are entangled (e > 0) for 
n > n cr it(d) and separable (e < 0) for n < n cr a(d). Then we change d and study the functional 
dependence between n cr n and d\ as d — ¥ oo we approach the continuum limit. Assuming that 
n c nt(d) is always finite (as it turns out, it is), we can expect three possible situations: 

1. — ^-j-jT — y oo for d — y oo. 

n crit (d) 

This would mean that, in the limit, regions of size L arbitrarily distant from each other can 
be entangled [remember that L is the fixed physical length, while D is given by (|28]l], 

— ¥ for d — ¥ oo. 



In this case, the physical distance D below which we can see entanglement would vanish for 
the given region size L, so no entanglement could be seen between separated regions. 

^J^y ->• C(L) for d — y oo, < C(L) < oo. 

In this last case, there exists a distance D(L) — C(L)L above which regions of size L are 
always separable, but below which they can be entangled if the appropriate profile is chosen. 



With the definition of the profiles that we use here, parameter D defined in the last section gives exactly the 
distance between the near ends of the regions. This choice is obviously the most suitable for the analysis. 
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As we are going to see, the data give strong numerical indication in favor of the third case. 

Let us consider in detail, as an illustrating example, the results for L = \[2 (the specific value 
is only chosen for numerical convenience). For fixed values of the separation between the blocks d 
and of the number of sites in each block n, we search for the optimal profile that maximizes the 
entanglement parameter (j4]). We do so for d fixed and increasing values of n, until we find e > 0. 
We repeat the procedure for d = 1, . . . , 16, so that finally we have, for each d, the smallest block size 
n cr it{d) such that the parameter e (maximized over all profiles) is positive. The inverse coupling 
constant scales with nasa _1 = l+^ [as required by (|25|) ]. The results are shown in Table [I] 



d = 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


ricrit(d) = 


2 


8 


14 


20 


25 


31 


37 


43 


48 


54 


60 


65 


71 


77 


83 


88 



TABLE I. Minimal block size n cr it that allows entanglement for fixed physical size of the regions L — v2 
as a function of the number d of oscillators separating the two regions. 



One can see that the relation between n cr i t and d is approximately linear, corroborating hypoth- 
esis 3 above. In Fig. [2] three optimal profiles are plotted, from which it can be clearly seen how the 
same shape is reproduced while increasing the number of points, as it is expected when approaching 
the continuum limit. We can conjecture that the continuous curve appearing in the limit n oo 
would correspond to the optimal profile in the continuum. 
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FIG. 2. Three optimal (maximizing the entanglement measure e) profiles for fixed physical size of the 
regions L — \f2. The number n cr it of oscillators in one block is the smallest that allows entanglement for 
the given number d of oscillators between the two blocks. For plot (a) d = 2, n cr u — 8, for (b) d = 8, 
n-crit = 43 and for (c) d = 16, n crit — 88. The plotted profiles are for the left block; the profiles for the 
block on the right have the same shape, mirrored. 



It is possible to extract the limiting value C(L) = lim^oo- — ^t^j by linear interpolation of the 

data in Table E We obtain C(L = \/2) = 0.17. From this, using relation (|28|) , we can calculate the 
physical distance below which regions of size \/2 can be entangled: D = C(L)L = 0.25 (expressed 
in units of the Compton wavelength). 

The same calculation can be repeated for different region sizes L. The linear behavior is confirmed 
for all the cases, as can be seen in the examples plotted in Fig. [31 

Once the linear dependence of n cr it on d is established, it is possible to calculate for each value of L 
the coefficient C(L). As can be seen from Fig. [H it is reasonable to infer that qjj\ has an asymptote 

for L — > oo, namely C(L) ~ in the figure the plot of (jjj^j is compared with a handmade 
asymptote with c — 1. This result implies that, as the region sizes become arbitrarily large, 
the distance above which they are separable converges to a finite value D cr u — lim£_ i . 00 Z?(L) = 
limi_j. 00 C(L)L ~ 1. The numerical results are stronger than our general proof from the preceding 
section as they imply the existence of a critical distance above which collective modes are separable 



14 



n Md) 
150 

100 



♦ 


#=1 


• 


y = 10 


■ 


y=50 



50 



• d 

5 10 15 

FIG. 3. Dependence of the critical number n cr it(d) of sites within the blocks on the number d of sites 
between the blocks for different region sizes L. The critical value n cr u(d) is defined as the minimal number 
of sites that gives entanglement of the corresponding collective operators. For each L a linear law is evident, 
the linear coefficient is defined as l/C(L). 



1 

coo 




L 



2 4 6 8 10 



FIG. 4. Dependence of the coefficient C(L) on the region size L. In the graph pjj-r is compared with a 
linear law. 



for any choice of the detection profiles, while the analytical approach assures only that for any 
pair of profiles there exists a critical distance for finding entanglement between the corresponding 
collective modes. It is remarkable that our calculation indicates that this distance is of the order 
of the Compton wavelength. 



V. RESULTS FOR THE MASSIVE KLEIN-GORDON FIELD 



Asymmetric triangular profiles. We will follow the intuition on the optimal profiles obtained 
from the linear harmonic chain. The most profound feature of the optimal functions found in that 
case is their asymmetry (see Fig. [2J and the fact that they are mirror images of one another: 
9a( x ) = 9b{~x). Thus, as a first approximation, to define the subsystems we use functions in 
the shape of asymmetric triangles. We restrict our numerical analysis to the one-dimensional case, 
however it is obvious that it may be extended to more dimensions. We parametrize each of the 
triangular profiles with the length of their support La/b and the position of the tip sa/b € (0, 1), 
where sa/b = \ gives a symmetric triangle. The property of mirroring results in both profiles 
having the same support size and tip position related by Sb = 1 — Sa, so that normalized profiles 
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are given by 

f/H'Sr) farxe(-L,-L(l-*)]; 

ffA(ac) = (s, L, x) := I ^.(^_^ for a: G (-£(1 - s),0); 9b(x) = g(s,L,-x) (29) 

[ otherwise. 

These functions not only satisfy condition (fT0|) but also assumption A3. The latter is evident once 
we realize that Fourier transforms of the profiles are related by complex conjugation, gA(k) = ffs(fc), 
which is a direct consequence of the fact that the triangles are mirror images of one another. Notice 
that the profiles are chosen in such a way, that parameter D = xb — xa is equal to the distance 
between their supports, D = meaning neighboring triangles. 



\ 0.8 




1 0.6 




\ 0.4 




\0.2 





-6-4-2 2 4 6 

FIG. 5. Two mirroring triangular profiles for the support's size L = 4, tip position s = 0.75 and distance 
between the supports D = 2. 



Moreover, again, we set field's mass m to 1 so that quantities of the length dimension are measured 
in Compton wavelengths. This finally makes the degree of entanglement depend only on three 
parameters: tip position s, size of the profiles' supports L and their distance D. An exemplary 
setting is shown in the Fig. [5l which in our scheme corresponds to the situation where two bosonic 
modes defined by the profiles are associated with the regions of length 4 separated by interval of 
length 2 in Compton wavelength units. 

With this choice of profiles, entanglement is found numerically for a certain range of parameters. 
Here we present a summary of our results. First of all, for separation D larger than D cr u 0.3 
the modes become separable - no entanglement can be found for any choice of the remaining two 
parameters. From our analysis of the linear harmonic chain, the critical distance was estimated to be 
of order 1. It was, however, done for optimal profiles to which triangles are just an approximation. 
For each separation D < D cr i t entanglement appears, once the size of the supports L exceeds 
some minimal value L m i n (D). This minimal length increases with D. The existence of the critical 
distance is manifested by L min (D) — > oo for D — > D crit ; see Fig. [51 

For given separation D, we can maximize the entanglement measure ^ over the size of the 
supports L and tip position s. In this way we obtain the maximal available entanglement as a 
function of the subsystems' separation D (Fig. [7]). For D — > D cr i t entanglement goes to zero. 
From the logarithmic plot we infer that in the intermediate range of separation parameter values, 
the entanglement measure e optimized over the two remaining parameters (L and s) decreases 
exponentially with the distance D. 

In Fig. [5] we show the dependence of the optimal (maximizing entanglement) values of the size 
of the supports L opt and tip position s op t- Qualitatively, the behavior of L opt and L m i n is the 
same. L opt is of the order of the Compton wavelength in the intermediate region of the separation 
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FIG. 6. Minimal length of the profiles' support Lmin, for which entanglement appears, plotted as a function 
of their separation D. Inset: logarithmic plot of the same dependence. 
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FIG. 7. (a) Degree of entanglement e maximized over the size of the profiles' support L and the tip 
position s as a function of the separation D, i.e. maximal available entanglement max[e] as a function of 
the separation of the profiles D. (b): Logarithmic plot of the same dependence. 



parameter values. More interesting is the dependence of the optimal tip position s opt on the 
separation D. It reaches its minimal value s op t ~ 0.84 for D « 0.2. Both for D — > and D — > D cr u 
the optimal tip position goes to 1. The difference between L opt and L m i n just confirms the result 



of [18j, that entanglement may emerge by going to larger blocks. 
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FIG. 8. Optimal and minimal length of the profiles' supports, plot (a), and optimal position of the triangles' 
tip, plot (b), as functions of the separation D of the profiles' support. By optimal values of parameters we 
understand such that maximize our entanglement measure e. 



For the sake of completeness we include the plots of entanglement as a function of the size of the 



17 



profiles' supports for separation D — 0.2 [Fig. 13(a)] and D = 0.1 [Fig.[9](b)]. To obtain each point 
of the plots, we optimized over the tip position. 
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FIG. 9. Entanglement as a function of the size of the profiles' support L for given separation of the profiles, 
D = 0.2 for plot (a) and D = 0.1 for plot (b). To obtain each point on the plots, the degree of entanglement 
e was maximized over the tip position s. 



The larger the values of the separation are and the farther the size of the profiles' supports is 
from its optimal value (for given separation), the more sensitive our entanglement measure is to 
changes of the tip position. In other words: as the separation D increases, e becomes more and 
more peaked over optimal values of L and s parameters. We exemplify this in Fig. 1101 where 
entanglement is plotted for separation D = [Fig. [10] (a)] and D = 0.15 [Fig. [10] (b)] as a function 
of the tip position for a few values of the size of the triangles' supports. 



(a)D "° (b)D=0.15 




FIG. 10. Entanglement as a function of the tip position s for a few chosen sizes of profiles' supports, L, 
and the separation of the profiles: D — for plot (a) and D — 0.15 for plot (b). L opt indicates the value of 
L parameter giving maximal entanglement for the relevant value of the separation parameter D, (compare 
Fig.©. 



At the beginning of this section, we stressed that the important feature of the profiles defining 
entangled modes is their asymmetry. Naturally, there arises the question to which extent this 
property is crucial. Is asymmetry necessary to obtain entanglement? In general the answer is 
negative, however, as shown in Fig. Qj] (a), the maximal amount of entanglement available in the 
case of symmetric profiles is more than 1 order of magnitude smaller than in general situation 
(i.e. when we can vary position of the tip of the triangles). Moreover, symmetrical modes become 
separable when the distance between their supports exceeds a value about 0.01 which is also more 
than an order of magnitude smaller than the critical distance estimated in the general, asymmetric, 
case (D cr it ~ 0.3). In Fig.QT](b) we plotted the dependence of the optimal size of symmetric profiles 
as a function their separation. There are also included points from the Fig. [5] (a), i.e. optimal lengths 
of the more general profiles for given separation. The optimal lengths of the symmetric profiles and 
of the general ones are, unlike the degree of entanglement, comparable. This result shows that the 
asymmetry of the profiles is indeed important for maximizing the degree of entanglement for given 
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distance between the subsystems. 
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FIG. 11. Results for the symmetrical triangles, i.e. profiles (|29|l with s = ^. Plot (a): for given separation D, 
we maximize e over the size of the supports L keeping s = | . In this way, we obtain the maximal available 
entanglement for symmetrical profiles as a function of the separation. Plot (b): optimal lengths of the 
symmetric and general profiles. By optimal we understand those that maximize the degree of entanglement 
e. 



We also stress that the exact shape of the triangle is not necessary. Since the space of compactly 
supported, smooth functions is dense in L P (R 3 ) for p > 1, all the integrals in the expression for 
e obtained with the triangles may be approximated to arbitrary precision with smooth functions 
satisfying all demanded constraints. Also an explicit "smoothing" is possible with the use of 

a series of Gaussian functions f a (x) = ^^ e_ ° 7 - They approach the Dirac delta in the weak 

limit i.e. lim a _>.o+ fa{x)g{x) dx = g(0). Convolving triangles with the function fa{x) for a 
sufficiently close to 0, we obtain a smooth profile that, again, approximates the entanglement 
measure e calculated with triangular profiles up to arbitrary precision. The price is that canonical 
commutation relations are also satisfied approximately - as Gaussians are not compactly supported, 
so are the final smoothened "triangles." 

Finally, we would like to point out that the joint state of the two regions, which we here consider, 
is in principle mixed, so it may no longer be useful e.g. for quantum communication. As proved 
in [33], for Gaussian two- mode systems entanglement (i.e. e > 0) is equivalent to distillability of 
the state. However, it is impossible to distill entanglement utilizing two identical, symmetrical 
copies of a two-mode Gaussian state at a time and performing Gaussianity preserving operations 
supported by classical communication [34j (non-Gaussian operations would be required). Entangle- 
ment distillation (or purification) is a transformation that consists of local operations and classical 
communication bringing several copies of a mixed entangled state into (approximately) pure entan- 
gled states which can further be utilized for quantum communication and quantum computation. 
(See Sec. IIIII for a justification that the state we define with the detection profiles is indeed Gaus- 
sian.) 



VI. CONCLUSION 



In this paper, we investigated entanglement between two regions of a scalar Klein-Gordon field 
in the vacuum state. By spatially integrating over field operators (and conjugate momenta) with 
two real, compactly supported functions we defined two pairs of collective phase space operators 
representing two bosonic modes, i.e. - subsystems corresponding to the regions given by the supports 
of the functions. Reducing the vacuum of a scalar field to these particular modes, we studied 
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entanglement between them with the use of an entanglement measure for two-mode continuous 
variables states (based on (l7|). 

For every two functions satisfying the aforementioned constraints (see Sec. it is proved 

that the corresponding subsystems are separable if the distance between them (i.e. between the 
supports of the functions) is larger than some finite value. From a numerical analysis of the 
discretized Klein-Gordon field we obtained strong indication that all considered modes become 
separable for separations larger than 1 Compton wavelength. We also gave an explicit example of a 
pair of functions (asymmetric triangles) that define entangled modes and investigated numerically 
the amount of entanglement in the corresponding system. 

The approach presented here aimed to take into account limitations that real experiments put 
on the properties that can be measured. From this point of view, our definition of observables 
is a reasonable first approximation towards a fully realistic treatment of the problem of vacuum 
entanglement as it assumes that only localized collective operators can be observed (since we cannot 
resolve field operators in single space-time points). Also an interaction has been proposed that could 
implement the desired measurements. 

According to our numerical analysis of the discretized field, in half spaces separated by more 
than a Compton wavelength we cannot find any entangled modes defined with the considered 
functions. In [l3 |. a violation of Bell's inequalities in the vacuum is in principle possible at arbitrary 
separations. Although our result cannot be directly compared with works on Bell's inequality 
violations (because our observables are Gaussian and on a Gaussian state such operators cannot 
show violation of Bell's inequalities [HI), we would like to point that the algebra of observables 
considered in 14] is much richer than ours. This implies that the operators needed for this violation 
(in the large separation regime) must be of a more intricate form than proposed here ( e.g. involving 
higher powers of the field operators). On the other hand, as the authors of [lH comment, if the 
distance D between the regions probed is much larger than a few Compton wavelengths of the 
lightest particle in the theory then the maximal Bell violation in the vacuum will necessarily be 
too small to be observed^- This shows that even with a wide range of allowed observables vacuum 
entanglement, if at all accessible, should be tested at small distances (of the order of the field's 
Compton wavelength^]). Our result asserts that, once we can probe close enough regions, to access 
vacuum entanglement we can restrict our measurements only to very simple observables, namely, 
field (and conjugate momentum) operators averaged over spatial regions. Such observables are 
usually considered as simplest to implement, at least in optics. 

There are several possible generalizations of the presented approach. First of all, each phase 
space operator {Qa/Bi Pa/b) may be defined via a different profile and the restrictions given by 
assumptions A1-A3 may be abandoned. Both these situations are not covered by our considerations. 
Further, we may allow for more general observables than field operators averaged over spatial 



regions, as e.g. in 14| or by considering space-time regions. Also investigating entanglement between 



excitations of local Hamiltonians (i.e. restricted to some space region) in the global vacuum state 
may be an interesting line of research. 

Our results have importance for investigating whether the vacuum of a quantum field has any 
operational meaning and if it may be accessed as an entanglement resource. Moreover, the method 
here presented is directly applicable to other than vacuum bounded energy states and also to 
systems described within nonrelativistic quantum field theory. The most extreme example of the 
latter (if the range of the correlations is considered), the BEC state, will be studied elsewhere. 



In their approach, the parameter describing the maximal violation of Bell's inequalities with operators from two 
local algebras, assigned to space-time regions separated by D, decays exponentially with mD (where m stands for 
the mass of the lightest particle in the theory). 

The fact that entanglement is only a short distance property of the vacuum is further confirmed by the fact that 
it is possible to transform the vacuum into a separable state by means of a nonlocal unitary involving only points 
at a Compton wavelength distance [36l| . 
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